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Abstract: The impact of Lorentz violation on the dynamics of a scalar field is investigated. 
In particular, we study the dynamics of a scalar field in the scalar- vector-tensor theory 
where the vector field is constrained to be unity and time like. By taking a generic form 
of the scalar field action, a generalized dynamical equation for the scalar-vector-tensor 
theory of gravity is obtained to describe the cosmological solutions. We present a class of 
exact solutions for an ordinary scalar field or phantom field corresponding to a power law 
coupling vector and the Hubble parameter. As the results, we find a constant equation 
of state in de Sitter space-time and power law expansion with the quadratic of coupling 
vector, while a dynamic equation of state is obtained for n > 2. Then, we consider the 
inflationary scenario based on the Lorentz violating scalar-vector-tensor theory of gravity 
with general power-law coupling vector and two typical potentials: inverse power-law and 
power-law potentials. In fact, both the coupling vector and the potential models affect the 
dynamics of the inflationary solutions. Finally, we use the dynamical system formalism 
to study the attractor behavior of a cosmological model containing a scalar field endowed 
with a quadratic coupling vector and a chaotic potential. 
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1. Introduction 

Scalar field theory has become the generic playground for building cosmological models 
related to particle physics, in particular for obtaining inflationary cosmologies which is one 
of the most reliable concepts to describe the early stage of the universe. The key property 
of the laws of physics that makes inflation possible is the existence of states of matter that 
have a high energy density which cannot be rapidly lowered. The Inflationary scenario Q 
relies on the potential energy of a scalar (inflaton) field to drive a period of early universe 
acceleration. It has been thought that the early universe could be well characterized by a 
series of phase transitions, in which topological defects could be formed Pj. In the context 
of the string theory, the natural values of the gauge and gravitational couplings in our 4d 
universe are explained by the dynamics of 'moduli' scalar fields [Q, ||]. Moreover, originating 
from the work of Sen H (see also ||), the possibility of the tachyon field being a candidate 
for the inflaton has been extensively studied. The tachyon action is of the Dirac Born 
Infeld form Q which leads to an equation of state interpolating between —1 at early times 
and at late-times. This suggests the possibility that the tachyon can play the role of 
the inflaton at early times and the dark matter at late-times. Serious difficulties, however, 
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plagues the tachyonic inflation j8). These include large density perturbations, problem with 
reheating and formation of caustics. 



Recent observational evidences especially from the Type la Supernovae ||, 1C] and 
WMAP satellite missions [|ll|], indicate that we live in a favored spatially flat universe 
consisting approximately of 30% dark matter and 70% dark energy. In the framework 
of the General Relativity this means that about two thirds of the total energy density 
of the Universe consists of dark energy, the still unknown component with a relativistic 
negative pressure p < —p/3. The simplest candidate for dark energy is the cosmological A- 
term. During the cosmological evolution the A-term component has the constant (Lorentz 
invariant) energy density p and pressure p = —p. However, it has got the famous and serious 
fine-tuning problem, while the also elusive dark matter candidate might be a lightest and 
neutral supersymmetry particle with only gravity interaction. For this reason the different 
forms of dynamically changing dark energy with an effective equation of state w < — 1/3 
were proposed, instead of the constant vacuum energy density. As a particular example 



of dark energy, the scalar field with a slow rolling potential (quintessence) [12] is often 
considered. The possible generalization of quintessence is a k-essence fl3"|| , the scalar field 
with a non-canonical Lagrangian. Any such behavior would have far-reaching implications 
for particle physics. However, recent theory of gravity with the Lorentz violation |l4|, |l5[ 
are proposed. 

More recently, authors in Ref. [|l(| explored the Lorentz violating scenario in the context 
of the scalar- vector-tensor theory. They showed that the Lorentz violating vector affects the 
dynamics of the inflationary model. One of the interesting features of this scenario, is that 
the exact Lorentz violating inflationary solutions are related to the absence of the inflaton 
potential. In this case, the inflation is completely associated with the Lorentz violation. 
Depending on the value of the coupling parameter, the three kind of exact solutions are 
found: the power law inflation, de Sitter inflation, and the super-inflation. 

The purpose of this paper is to study the dynamics of a scalar field in the framework 
of Lorentz violating scalar-vector-tensor model, taking into account the effect of the dy- 
namically coupling vector. In this framework, we explore a class of exact solutions such 
as evolution of a scalar field and equation of state parameters. We discuss an inflationary 
scenario with a power-law coupling vector model with two typical potentials: an inverse 
power-law potential and a power-law potential. Then, we show that it is possible to find 
attractor solutions in the Lorentz violating scalar-vector-tensor model in which both the 
coupling function and the potential function are specified. 

The organization of this paper is the following: in Section ^, we set down the general 
formalism for the scalar-vector-tensor theory where the Lorentz symmetry is spontaneously 
broken due to the unit-norm vector field. In Section [3], we use our formalism to find an 
exact solution of the equation of state. In Section |], we study Lorentz violating (inverse) 
power law inflation. In Section ||, the critical points of the global system and their stability 
are presented. The final Section is devoted to the conclusions. 
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2. General Formalism 



In the present section, we develop the general reconstruction scheme for the scalar-vector- 
tensor theory. We will consider the properties of general four-dimensional universe, i.e. the 
universe where the four-dimensional space-time is allowed to contain any non-gravitational 
degree of freedom in the framework of Lorentz violating scalar-tensor-vector theory of 
gravity. Let us assume that the expectation values of a vector field is < Oju^u^lO >= —1. 
The action can be written as the sum of three distinct parts: 

S = Sg + Su + S,/, , (2.1) 

where the actions for the tensor field S g , the vector field S u , and the scalar field S"^ are, 
S 9 = jd^ 9T ^R, (2.2) 

S u = J d 4 xy/=g~ [-ft - (3 2 V^u u V vUfl - 3 (V M u'<) 2 

-ft//V/V A + A(u\|l)] , (2.3) 
S,p = I d 4 x^9 ^ ■ (2.4) 



In the above (3i(<j>) (i = 1,2,3,4) are arbitrary parameters which has the dimension of 
mass squared. It means that \J~f3l gives the mass scale of symmetry breakdown. is the 
Lagrangian density for scalar field, expressed as a function of the metric g^ u and the scalar 
field <j). Then, the action ( p.l| ) describes the scalar-tensor-vector theory of gravity. The 
dimensionless vector field, u^, satisfies the constraint 

u^Uf, = -1. (2.5) 

For the background solutions, we use the homogeneity and isotropy of the universe 
spacetime 

ds 2 = -M 2 {t)dt 2 + e 2a( - t) 5 ij dx i dx j , (2.6) 

where N is a lapse function and the scale of the universe is determined by a. We take the 
constraint 

1 



xr , 0,0, 1 , (2.7) 

where M = 1 is taken into account after the variation. Varying the action ( |2.1| ) with respect 
to g^ v ' , we have field equations 

iV - \d^R = 87rGr^ , (2.8) 

where = T$ + t]$, is the total energy- momentum tensor, T^/ and T^P are the 
energy-momentum tensors of vector and scalar fields, respectively, defined by the usual 
formulae 

T $ = "2^7 + 9^ k \ k = u,<f>. (2.9) 
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The time and space components of the total energy-momentum tensor are given by 



To = -p u -p<f> , T- = p u + ity , 
where the energy density and pressure of the vector field are given by 

-3fiH 2 , 



Pu 
Pu 



3 + 2— + 2^) f3H 2 



/5 = /?i + 3/3 2 +/3 3 • 



(2.10) 

(2.11) 
(2.12) 
(2.13) 



From the above equations, one can see that (5^ does not contribute to the background 
dynamics. A prime denotes the derivative of any quantities X with respect to a. X' 
is then related to its derivative with respect to t by X' = (dX/dt)H~ l = XH^ 1 where 
H = da I dt = a is the Hubble parameter. From Eqs. ( [2.11D and (|2.12 ), one obtains the 
energy equation for the vector field u 



Pu + HPu+Pu) =+3#V , 



(2.14) 



and for the scalar field 



(2.15) 



The total energy equation in the presence of both the vector and the scalar fields is, 
accordingly, 

p' + 3{p + p) = , (p = p u + P<j> ). (2.16) 

This energy conservation equation can also be obtained by equating the covariant diver- 
gence of the total energy-momentum tensor to zero, since the covariant divergence of the 
Einstein tensor is zero by its geometric construction. It follows from contraction of the 
geometric Bianchi identity. 

Substituting Eq. ( 2.1 0|) into the Einstein equations (E^§1 ), we obtain two independent 
equations, called the Friedmann equations, as follows: 



-3H 2 = 8irG (3f3H 2 



-2HH' - 3H 2 



8irG 



H' (5' 
3 + 2 lT + 2 J 



These Friedmann equations can be rewritten as 



1 + 
1 + 



1 



1 



8irG(3 J H W H 



(HH' + H 2 ) 



8ttG(3 



1 (P± 3 P<P 
"6 V/3 



9 



3 



(2.17) 
(2.18) 

(2.19) 
(2.20) 



The second term on RHS of Eq. ( 2. 20]) is a consequence of the coupling vector field as a 
function of scalar field. If fa = 0, thus without the vector field, the above equations reduce 
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to the conventional ones. And in the case (3 = const., the above equations are lead to the 
Friedmann equations given in Ref. [17]. 

Using Eqs. ( gig ) and (|2l5|) , we obtain a set of equations as follows: 



H' & 3 , 



, 



^-^-2 (1+ ^ ) = °' 



where 



and oj^ 



^ + f + 3(l+^ 



= P + 







(2.21) 
(2.22) 
(2.23) 

(2.24) 



is the equation of state of the scalar field. It is easy to check that the 



equations ( 2.21 )-( 2.23 ) satisfy the following constraint 



H P4> 







(2.25) 



In order to solve the Eqs. (2.21)-(2.23) and ( 2.251 ), we have to specify the model and the 
matter content of the universe. The general solution of these equations can be written as 



H(3 oc exp 
H 

— oc exp 

P<t> 

oc exp 



-(1 + u^{a))da 



-(1 + LO <f) (a))da 



3(1 + u}^(a))d> 



a 



(2.26) 
(2.27) 
(2.28) 



If the functions ujs and /3 are given, then we can find the evolution of the Hubble parameter 
under the Lorentz violation. For example, the cosmological constant corresponds to a fluid 
with a constant equation of state uj^ = —1. Thus the above equations reduce to: H(3 oc 1, 
H cc and p^fi oc 1 where H, p^ and j3 are functions of a. If oj^ is a constant parameter 
of a simple one component fluid, and for a given a(t), Eqs. ( 2.2l|) -( 2.23| ) can be used to 
determine (3(a) and p^(a). We, then, are able to determine the potential of the Lorentz 
violation model. 



3. Dynamical Equations for Scalar Fields 

For a given scalar field Lagrangian with the FRW background, we can obtain the equations 
of motion for a scalar field by using Eq. ( g3|) and Eqs. Q2.2Ij )- (^2^) . Let us consider the 
Lagrangian density of a scalar field <§> with a potential V((j)) in Eq. ([2.1|): 



(3.1) 
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where (V(/>) 2 = g^d^d^cp. Ordinary scalar fields correspond to r] = 1 while rj = — 1 is for 
phantoms. For the homogeneous field the density p^ and pressure p^ of the scalar field, 
may be found as follows 



P<f> 



|tfV 2 + V{4>) , 

_ v{(j)) 



The corresponding equation of state parameter is, accordingly 

= 1 - r]H 2 <P' 2 /2V 

^ P4, 1 + V H 2 <P> 2 /2V ' 

Substituting Eq. (|3.2|) into Eq. (|2.19|) , the Friedmann equation leads to 



3/3 



1 r^2 / 2 + y(0) 



(3.2) 
(3.3) 

(3.4) 
(3.5) 



Now, differentiating Eq. ( |3.2| ) with respect to a and using Eq. (|2.15| ), and also differentiating 
Eq. ( |3.5| ) with respect to a and using Eq. ( |3.6|) give, respectively, 



If 



+ 3 U' - rj 



Ya 

H 2 



(3.6) 
(3.7) 



Substituting Eq. 
written as 



into the Friedmann equation the potential of the scalar field can be 



V = 3/3H 2 



3^ It 



+ 



H,4 

IT 



(3.8) 



Note that in the above equations the Hubble parameter H has been expressed as a function 
of <f>, H = H((j){t)). From Eq. (2.21), the equation of state can be written as 

J , / «... (.,, V 



i 



, 0' 2 

-1 + -77_ 

3 ' p 



(3.9) 



Equations ( |3.7| ) and ( |3.9[ ) are two equations that we need to solve for the scalar field <$> 
and the equation of state wq. This is achieved only if the Hubble parameter H{(j>) and the 
coupling vector (3(4>) are known. For different choice of the Hubble parameter H((f>) and 
the coupling vector /3(0), it is possible to extract a class of exact solutions of Eqs. Q3.7p 



and ( |3.9| ). We shall solve Eqs. (3/7) and ( |3.9| ) to obtain the following physical quantities 
(V and K are the potential and kinetic energies, respectively): 



V = -(l-^)PH 2 



p^ = 3f3H , p^ = Suj^PH . 

In the following two subsections we will explore a class of exact solutions. 



(3.10) 
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3.1 Exact solutions and the behavior of scalar fields 



We shall have to solve equations ( |3,7p and (3.9) for H, uj^, (3, and V, which is not possible 



unless two are known. In the present subsection, we consider an example to find an exact 
solution of the equation of state of the scalar field in the quadratic coupling vector. The 
equation of state for the scalar field has been intensively studied in |l8| for the so called 
tracking cosmological solutions introduced in and some classes of potentials allowing 
for the field equation of state were described. 
Let us consider a simple model 



H = H , 0(4) = mtf , (3.11) 

where Hq and m are positive constant parameters. The equation ([^7]) can now be inte- 
grated to yield the evolution of the scalar field 

<f>(t) = <f> exp[-4r l rnH (t-to)] , (3.12) 

where <p(t = to) = <fio is a constant. Then, it is easy to find the equation of state of the 
scalar field by using Eq. (^S|). We obtain 



X6 

uj^ = —1 + i f or ordinary scalar fields , (3.13) 

X6 

ujcf, = — 1 — — m , for phantom fields . (3-14) 

Then, the potential and the kinetic energies, the energy density and the pressure of the 
scalar field evolve according to 

V(t) = mH^l(3 - 8m) exp [-8rjmH (t - t )] , (3.15) 

K(t) = 8r/ (mi/ o 0o) 2 exp [-8r}mH (t - to)] , (3.16) 

p(t) = 3mH^exp[-8r)mH (t-t )] , (3.17) 

p(t) = r]mHl<t>l (16m - 3r?) exp [-8r]mH (t - t )] . (3.18) 

The above solutions are completely associated with the Lorentz violation. The model 



( 3.11 ) depicts that the cosmic evolution starts from a constant value of the scale factor and 
grows exponentially, a(t) = aoe H °^~ to \ The coupling vector decreases exponentially for 
the ordinary scalar field and increases exponentially for the phantom field from a constant 
value of m</>Q. Hence the potential energy, the kinetic energy, the energy density and the 
pressure decrease exponentially for the ordinary scalar field. For the phantom field, on the 
other hand, the potential energy and energy density as well as the absolute values of the 
kinetic energy and the pressure increase exponentially. Note that the kinetic energy and 
the pressure begin with the negative values. The Eqs. (3.13)-(3.14) show that the equation 



of state u)<f) is non-dynamical because it only depends on the value of the coupling vector 
parameter m, both for the ordinary scalar and the phantom fields. Since an accelerated 
expansion occurs for tu^ < —1/3 then we have m < 1/8 for the ordinary scalar field. 
However, the present data of the Universe seems to tell that oj^ might be less than — 1. 
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Thus, the value m may be chosen in order to fit the present observable constraint on the 
equation of state parameter. 

In other case, for instance, H(4>) = Hq^ and j3(4>) = mcf) 2 , we also find the constant 
equation of state, 



city = -1 + -r/m(£ + 2) 2 . 

The condition for the accelerating Universe d or H' / H > — 1 yields 

1 



This model gives a power law expansion 



2£(£ + 2) 



(3.19) 



(3.20) 



a(t) 



i + _2ro ( 



+ —^{t-t ) 

p 



p > 1 



where 



The scalar field evolve as 



1 



P 



2r)m£(£ + 2) 



(3.21) 



(3.22) 



(3.23) 



Hence, the complete set of solutions is found by substituting Eqs. (3.19) and (3.23) into 
Eqs. ( pCl) . 

In the following subsection, we will see that the equation of state may be dynamics. 
For this purpose we generalize the coupling vector to j3((f>) = m<j) n , n > 2. 

3.2 Variable equation of states 

Let us consider a model where the coupling vector is a power law of the scalar field, 



H = H , p{<j)) = mtf 1 , n > 2 , 



(3.24) 



where Hq, m and n are constant positive parameters. Following the same above procedure, 
the scalar cfi can be evaluated as, 



<f>(t) 



[1 + 2rimnH (n - 2)#J~ 2 (t - t Q ) 
the coupling vector is given by 



1 ! 

n— 2 



(3.25) 



[1 + 2rimnH (n - 2)^~ 2 (t - t )] 



n i 
n — 2 



(3.26) 
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and the dynamical equation of state (^^) is 

uj(t) = -1 + -2 — '—^ . (3.27) 

l + 2T/mnflb(n-2)^- 2 (t-to) 

Then, the potential and kinetic energies, the energy density and the pressure of the scalar 
field are given by 

2r]mn 2 (j)^ 2 /3 



V(t) = 3mH { 



1 



1 + 2-qmnH (n - 2)(^' 2 {t - t ) 
1 



[1 + 2r)mnH {n - 2)^- 2 {t - to)] ' 



(3.28) 



K(t) _ 2 "W' 2 ^ , (3.29) 

[1 + 2rjmnH (n -2)<f%- 2 (t- t )} n ~ 2 

p(t) = ^ — , (3.30) 

[1 + 2r]mnH (n - 2)<^" 2 (t - t )] - 2 

4r/mn 2 ^~ 2 /3 



p(t) = 3mH 2 ^ 



-1 + 



1 + 2r]mnH (n - 2)^~ 2 (t - t ) 
1 



. (3.31) 

[1 + 2r]mnH (n - 2)^~ 2 (t - t )] ™" 2 

Thus, the model fl3.24|) describes that the cosmic evolution grows exponentially from 
a constant value of the scale factor, a(t) = aoe H °( t ~ t °\ while the coupling vector /3 started 
from a constant value of the scalar field, m(j>Q. The equation of state is dynamical both 
for the ordinary scalar and phantom fields. Then the potential energy, kinetic energy, the 
energy density and the pressure decrease for the ordinary scalar field. For the phantom 
field, on the other hand, the potential and energy density increase while the kinetic energy 
and pressure begin with the negative values. 

4. Lorentz Violating Inflation Scenario 



As it has been studied by authors in Ref. [16|, the Lorentz violation on the inflationary 



scenario can be divided into two parts: the Lorentz violations stage SttG(5 3> 1 and the 



standard slow roll stage 8irGf3 <C 1. The first stage corresponds to /3 = (3 in Eq. ( 2.24 ) and 
the second stage corresponds to (3 = l/8nG, then we have the usual dynamical equations. 

In this section we will consider the inflationary scenario for the scalar field (inflaton). 
In particular, we consider a power-law coupling vector, 0(<j>) = m<j) n , with two types of 
the potential: V(<p) = /i 4+1/ _ly and V((f>) = ^M 2 (f> 2 . Here [A, v and M are parameters. 
Thus, the dynamics of each particular inflationary model are determined by the Friedmann 
equation and the scalar field equation of motion once the functional form of the inflaton 
potential and the coupling parameter have been specified. Let us collect the dynamics- 
related equations for the inflaton the Friedmann equation (|3.5| ) in inflationary models 



H 



2 



1 

3/3 



l -H 2 ^' 2 + V(4>) 



(4.1) 
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the constraint equation (obtained from Eqs. (]3.4|) and ( 2.21|) ) 



^ + i£ + £ = 0, (4.2) 



and the equation of motion Eq. 

4>" + ^ + 30' + ^ + 3^ = . (4.3) 

/3 is given by Eq. (|2.24j) . Then, at the critical value of cf>, the effective coupling vector 
becomes 

8irG(3(<f> c ) = 1 . (4.4) 
For example, a coupling parameter of the form f3 = mcf) 2 gives the critical value 

Mpl M v i = G^. (4.5) 



ran 



Let (j>i be the corresponding initial value of the scalar field. Putting <pi ~ 3M p i, the Lorentz 
violation implies the criterion m > l/(727r) ~ 1/226. 



The set of Eqs. (4.1)-( |4.3| ) constitutes the equations we have to solve for the problem 
specified by the coupling parameter /3(<f>) and the potential V(cfi). In following subsection, 
we consider with a model with the coupling parameter f3{4>) is given by 

(3(cP) = mr , (4.6) 

where n and m are parameters. For the model (4.6), we obtain the critical value of the 
scalar field and the criterion for Lorentz violation 

*■-(*£) mi m> ^w- (47) 

Now, we consider two typical potentials appear in many cosmological implications: an 
inverse powerdaw potential and a powerdaw potential. We discuss those solutions and 
analyze the two regimes separately. 

4.1 Inverse power law potential: V(<p) = fi +v (j)~ L ' 

In this subsection, we consider the class of power law potential 

V{<}>) = //+>-" , (4.8) 

where and v are constants. Inverse power law models are interesting for a number 
of reasons. In conventional cosmology, they drive intermediate inflation [20] and typically 
produce significant tensor perturbations for almost scale-invariant scalar fluctuations. They 
arise in supersymmetric condensate models of QCD [21] and can in principle act as a source 
of quintessence [12, 22 1. 
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4.1.1 Lorentz violating stage 

Let us first consider the Lorentz violating stage, 8irG(3 3> 1 {J3 = ft) , we have the equations 
(|4.1j) — Q4.3|) . In this stage both the coupling function and the potential function are relevant. 

An inflationary epoch, in which the scale factors a are accelerating, requires the scalar 
field <j> to evolve slowly compared to the expansion of the universe. Thus, the following 
conditions of slow-rolling are required: 



H 2 6' 2 < V 



" < (j)' , </>' 2 < ft, and ft'<ft 



(4.9) 



The formalism which gives these slow roll conditions are discussed in Ref. (jOJ. This is 
sufficient to guarantee inflation. Under the slow-roll conditions Eq. fl4.9|) , the Eqs. (|4.l| )- 
(|4.3|) can be simplified. We obtain the slow roll equations 



V 



P,4> , V ,<t> 



H 2 ~ ^ , and </>' ~ -ft + ^ \ . (4.10) 
Inserting Eq. ( f4.6[ ) and the potential of the form V(<p) = ix i+u (f)~ u into Eq. (4.10), we have 



h 2 = v — r (u+n) , 

3m 

^ = - m (n - v)^ n "^ 



One can then solve for <p from Eq. fl4.12|) , 

4>{a) = [4>1~ n + m(n - 2)(n - v){a - an)] 



i 

' n-2 



for n ^ 2, n ^ v 



where d>(a = at* 



= fa is a constant. The Friedmann equation gives 

# («) = V~ [0i" n + m(n - 2)(n - i/)(q - a;)] - 2 . 
6m 



(4.11) 
(4.12) 

(4.13) 
(4.14) 



The solution ( [4.13 ) and the slow roll conditions ( [4.9|) during the Lorentz violating stage 
give n > 2 because 



,/2 I ^ a -2(l-n)/(2-n) 



ft 



-2(l-n)/(2-n) 



«.ft 



a n/(2-n) 
n/(2-n) 



(4.15) 
(4.16) 



From Eq. (4.14), the universe expands during the Lorentz violating stage as 



o(t) 



exp 



5 1 

+ 



1 



c c v^- 1 

where the constants A,B,C and D are 



AC(-D-l)(i-tj 



i 

' D-l 



, c/o 



(4.17) 



A 



3m 



B 



k 2-n 



, C = m(n — 2)(n — z/) 



D 



n + v 
2{n-2) 



(4.18) 
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Combining Eqs. ( 4.17 ), ( 4.13| ) and ( [4.14 ), we obtain the physical quantities 



<P(t) = [ W Zi-AC(D-l){t-t l ] 



(£>-l)(n-2) 



D 

'(D-l) 



H(t) = A[ WJ - I - AC(D -l)(t-U] 



1 \ (D-l)(n-2) 

m=n\ wz - 1 -AC(D-l)(t-t l ) 



(4.19) 
(4.20) 
(4.21) 



The scalar field energy density, on the other hand, evolves according to 

1 



p(t) ~ V = 3mA^ 



B D 



_-AC(D-l)(t-U] 



n-2D(n-2) 
(D-l)(n-2) 



(4.22) 



One can see that the Hubble parameter H decreases during the Lorentz violation stage. 
For n = 2, v ^ 2, Eq. ( |4.6j ) and the second part of Eq. ( [4.10 ) gives 



(4.23) 



where 0(a = aj) = 0j. For this solution to satisfy slow roll conditions ( |4.9| ), we need 
m < 1/(2 — z/) 2 . Thus, we have the range 1/226 < m < 1/(2 — v) 2 of the parameter for 
which the Lorentz violating inflation is relevant. The Hubble parameter as a function of 
the scale factor, a, is given by 



H 2 (a) = ^2 e -m(^-4)(a- ai ) 



where 



H 



4+!/ 



3m(f>. 



(2+u) 



(4.24) 



(4.25) 



while the scale factor a(t) = e a is of the form: 



a(t) 



0(t) 



i(i/-2) 



(4.26) 



Now we obtain the evolution of some physical quantities as follows 

a(t) 

0(t) 



rruf) 2 
P(t) = 



[m(u 2 -A)H 2 {t-ti)]™^-*) , 


(4.27) 


[m(i/ 2 -4)ff?(t-ii)]^ ) 


(4.28) 


= ^m{A-v 2 ){t-ti) , 


(4.29) 


= [m{u 2 -A)H 2 (t-U)]^ , 


(4.30) 


^ [ m (i, 2 -4)# 2 (i -£;)]" ^ • 


(4.31) 
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4.1.2 Standard slow roll stage 



The governing equations (4.1)-( |4~3|) in the standard slow roll stage 87rG/3 <C 1 (j3 
(87rG) _1 ), are, accordingly, 



8itG 



-H 2 6' 2 + V 



^ + 4vrG0' 2 = O, 

0" + §0' + 3*' + ^ = O. 



In this case the slow roll equations are given by 



8vrG V 



(4.32) 

(4.33) 
(4.34) 

(4.35) 



For the potential model V(4>) = fi A+u (j) 1 ' i the evolution of the inflaton and the Hubble 
parameter can be solved as 



2 (a) = (f? c + - — —(a — a c ) 



4vrG 



if 2 (a) 

and the scale factor is given by 

a(t) 



8ttG 



4+u 



+ 



4vrG 



(a - a c ) 



-u/2 



(4.36) 
(4.37) 



exp 



4vrG 



(0 2 (t)-<A 



(4.38) 



The evolution equations are given by 



— = ex P {-P- + ?r + A S C S {D S + l)(t - t c )] ^ 

<P(t) = [B^ +1 + A S C S (D S + l)(t - t c )] ^^tt j 
H(t) = A s [B°° +1 + A S C S (D S + l)(t - t c )] ^ , 



and the scalar field energy density evolves as 

2D S 

[B°° +l + ^ S C S (D S + l)(t - t c )] 



3 ^ 2 C S 



8 V A 

where ^4 S , B s , C s and D s are the constants 



(4.39) 

(4.40) 
(4.41) 

(4.42) 



A* 



8itG 



3 r- . rc > 47rG . 4 

Note that, in the standard slow roll stage, the Hubble parameter H increases. 



(4.43) 
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Another interesting quantity is the number of e-folding during the inflationary phase. 
The total e-folding number reads 

B 1/1 

N = ~C + C " AC{D ~ l){tc ~ k 

-T^ + TT [B? s+1 + A S C S {D S + l)(t e - t c )] 
= £ (^" n - $T n ) + ^ " <t>l) , (4-44) 



for m > 2, v ^ m and 



N = n{u 2 _ 4) lQ g Hv 2 -4)H?(t c - tl )] 

-T^ + TT l B ? S+l + A sCs(Ds + l)(te ~ Q] 

= -7^—, log (I 1 ) + — W - $ ■ ( 4 - 45 ) 

77i(2 — Z/j \<Pc/ ^ 

for n = 2, z/ ^ 2. Note that the first terms of the above equations arise from the Lorentz 
violating stage. As an example, let us take the values: N = 70, m = 10 -2 , n = 2 and 
v = 1. If cf) e ~ 0.3Af p ; is the value of scalar field at the end of inflation, then, <p c ~ 2M p i. 
The contribution from the inflation end is still relevant. Therefore, we get <pi ~ 2.5M p [. 

4.2 Power law potential: V{<f>) = \M 2 4> 2 
4.2.1 Lorentz violating stage 

The most realistic inflationary universe scenarios are chaotic models. For the model V{<j)) = 
^M 2 (f) 2 , assuming the slow roll conditions, we find the slow roll equations during the Lorentz 
violating regime as follows 

#2 = Ml r (n-2) (446) 
6n 

(/)' = - m (n + 2)<^ n ~ 1 ) . (4.47) 
Then we find the solution (|4.47|) as 

0(a) = [cp 2 ~ n + m{n 2 - 4) (a - a;)] ^ , (4.48) 

for m ^ 2 and 

^(a) = ^e- 4m(Q -" !) , (4.49) 



for n = 2. The inflationary scenario of this model was already obtained in Ref. |16[ 
where the Hubble parameter becomes constant during the Lorentz violating regime and 
1/226 < m < 1/16 is the range of parameter m. We concern here the solution for 
The solution for the Hubble parameter is given by 

M 2 

H 2 ( a ) = — [c/ ) 2 ~ n + m(n 2 -A)( a - ai )] , (4.50) 
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and 



a(t) 



exp 



m(n 2 - 4) \4> n - 2 (t) $ 



n-2 



(4.51) 



which is the solution for the scale factor. As in the previous subsection, we also obtain 
n > 2 which the effect of Lorentz violation occurs in this regime. The time evolution of 
the above equations can be obtained by integrating Eq. ( 4.5C| ), we get 

h i r . ,„ i 1 2 

a(t) = a.i 



b 1 

H 1 



b 1/2 + -dc(t - U 



where 



b = $~ n , c = m(n 2 - 4) 
Then the evolution equations are given by 



M 2 
6m 



n > 2 . 



a(t) 



exp 



b 1 

- + - 

c c 



b 1/2 + -dc(t - ti 



H(t) = d 
f3{t) = n 
p{t) = 3nd 2 



b 1 ' 2 + -dc(t - U) 



2 

' n-2 



b 1 ' 2 + ldc(t - U) 
b 1 ' 2 + ^dc(t - ti) 



b 1 ' 2 + -dc(t - U) 



2n 
' n-2 



4 

" n-2 



(4.52) 
(4.53) 

(4.54) 

(4.55) 
(4.56) 

(4.57) 

(4.58) 



Since b, c and d are positive constants, one can see that the Hubble parameter H and the 
scale factor a increase during the Lorentz violating stage for n > 2. In the case n = 2, 
the Hubble parameter is constant. In the following subsection, we will see that the Hubble 
parameter decreases in the standard slow roll stage. 

4.2.2 Standard slow roll stage 

Now, let us consider the chaotic inflationary scenario in the standard slow roll stage. A set 
of the dynamical equations of the scalar field are given by Eqs. ( [4.32| )-( |4.34| ). Assuming 
the standard slow roll conditions, we find the slow roll equations 



H 2 ~ ^-M 2 6 2 



1 



AttG 

The evolution of the inflaton can be solved as 

1 

2^G 



(a - ct c ) 



(4.59) 
(4.60) 

(4.61) 
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The Hubble parameter and the scale factor a(t) = e a can be also obtained as 



AttGM 2 



1 



(q - a c ) 



3 V TC 2 ^G 
a(t) = a c exp [27rG(<#? - (j) 2 (t))] 



From Eq. ( |4.62| ), we obtain 



b s 1 

a[t) — a c = 

C« Co 



b 1 / 2 - -d s c s (t-t c 



1 2 



and the dynamical evolutions are given by 



a(t) 



exp 



c.« c 



- ( bV 2 - ^d s c s (t - t c ) 



<t>{t) = b 1 / 2 - ^d s c s (t - t c ) , 
b l J 2 - \d s c s {t-t c 



H{t) = d 

3c s d, 



P(t) 



b 1 / 2 - -d s c s (t-t c ) 



1 2 



with 



1 



4 



AttGM 2 



a " c ' " 2vrG ' ~ a 3 
Note that the Hubble parameter decreases in the standard slow roll stage. 
In the case of chaotic potential, the total e-folding number reads 

1 2 



N 



b 1 
— + - 

c c 

1 



b 1 ' 2 + ^dc{t c - ti 



_ls _ 1_ 

Co Co 



b 2 ~ n - f-" 1 



) +2ttG 



(4.62) 
(4.63) 

(4.64) 

(4.65) 

(4.66) 
(4.67) 

(4.68) 
(4.69) 



(4.70) 



m(n 2 — 4) 

where (fi e is the value of scalar field at the end of inflation. Notice that the first term arises 
from the Lorentz violating stage. 



5. Phase-space analysis 

In this section, we investigate the global structure of the dynamical system via phase plane 
analysis and compute the cosmological evolution by numerical analysis. Introducing the 
following variables: 



V 



Ai 



.it 

PP,cf><t> 

PI ' 



y 



3H 2 (3 



Ya 

V 



r 2 = 



vv. rt 

V 



2 ' 



(5.1) 
(5.2) 
(5.3) 
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the Eqs. (|4.2| ) and ( [4.3|) can be written as a plane-autonomous system 



-3x(l - x 2 ) + A / -(Ax + A 2 )y 2 , 



3x 



(Ai + A 2 ) 



a; = -VEx 2 (r x - 1) 



a; 



-v^A 2 , 



2A 2 



(5.4) 
(5.5) 

(5.6) 
(5.7) 



where the prime denotes a derivative with respect to the logarithm of the scale factor, 
a = In a. The functions X\((j)) and A 2 (<^>) determine a type of the coupling vector and the 



potential, respectively. The Friedmann constraint, Eq. (4.1), takes the simple form 



x 2 + y 2 = 1 



(5.8) 



The equation of state for the scalar field could be expressed in terms of the new variables 
as 



W0 



y 



x 2 + y 2 



(5.9) 



Notice that x 2 measures the contribution to the expansion due to the scalar field kinetic 
energy and the coupling function, while y 2 measures the contribution to the expansion due 
to the potential energy and the the coupling function. 

Equations (5.4)— (5.7) are written as an autonomous phase system of the form x' = f(x) 
where x = (x, y, X\, A 2 ). The use of this form for the dynamical equations allows the fixed 
points of the system to be readily identified, and the so-called critical points xo are solutions 
of the system of equations f(xo) = 0. To determine their stability we need to perform 
linear perturbations around the critical points in the form x = xo + u, which results in the 
following equations of motion u' = Mu, where 



Mi. 



dxj 



(5.10) 



In the case of the dynamical equations (|5.4j )-( |5.7D , u is a 4-column vector consisting of 
the perturbations of x, y, X\ and A 2 . Thus, Mij is a 4 x 4 matrix. The stability of the 
critical points is determined by the eigenvalues /Uj of the matrix M at the critical points. A 
non-trivial critical point is called stable (unstable) whenever the eigenvalues of M are such 
that Re(fii) < (Re(fii) > 0). If neither of the aforementioned cases are accomplished, the 
critical point is called a saddle point. 

In the following, we will study the simplest model, 



(3(4>) = mcj) 2 



V(<t>) 



1 



-M 



2^,2 



(5.11) 
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Figure 1: The phase plane of Lorentz violating kinetic dominated solution for m > 3/8. 



where m and M are parameters. Substituting Eqs. ( 5.11 ) into Eqs. ( |5.2[) and (|5.3| ), respec- 
tively, we obtain 



Ai — Ao 



-2Jm 



Ti = r 2 = - 
i i 2 



(5.12) 



and Eqs. ( |5.6|) and (|5.?1) are trivially satisfied. In the former, Eqs. ( |5.4j) and (|5.5|) can be 
fused into the single equation, 



3x 



2^1 + A 2 ) 



(1-x 2 ) 



3x + 2v 6m ) (1 — x 



(5.13) 



which is one dimensional phase-space corresponding to the unit circle. Critical points cor- 
respond to fixed points where x' = 0, and there are Lorentz violation self-similar solutions 
with 



H' 

— = —3x 2 + VQX\x 
H 



(5.14) 



Note that the second term arises from Lorentz violation. Applying the above proce- 
dure, setting x' = 0, the critical points (xo,yo) of the system are (1,0), (—1,0), and 
{-^Sm/2>, - 8m/3). For any form of the potential in the Lorentz violation stage, the 
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critical points (1,0) or (—1,0) correspond to two Lorentz violation kinetic-dominated solu- 
tions. Then, the critical point (— y8m/3, yl — 8m/ 3) corresponds to a Lorentz violation 
potential-kinetic solution. Integration of Eq. (|5.14| ) with respect to a will show that all 
critical points, xq, correspond to the Hubble parameter 

Focexp^-^ . (5.15) 

This relates to an expanding universe with a scale factor a(t) given by a{t) ~ t p , where 




Figure 2: The phase plane of Lorentz violating kinetic-potential solution for m < 3/8. 



3xq — \/6\\xo 3xq + 2\/6mxo 

The linear perturbation about the points xo+ = +1 and xq- = — 1 give the eigenvalues 
//+ = 6 + A\/6m and /i_ = 6 — 4\/6m, respectively. Thus for positive m, xq + = +1 is 
always unstable and xq- = —1 is stable for m > 3/8 but unstable for m < 3/8. Moreover, 
in the linear perturbation about the Lorentz violation potential-kinetic solution, we obtain 
the eigenvalue // = 8m — 3. The solution is stable for m < 3/8. In Figs. |l] and |2[ we show 
the phase plane plot for m > 3/8 and m < 3/8. We note that the trajectories are confined 
inside the circle given by x 2 + y 2 = 1. 
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Another remarkable feature of the above model is that the equation of state is given 

by 

16 

uty = -1 + — m , (5-17) 

completely determined by the parameter m of the coupling vector. Thus, we always have 
> — 1 for ordinary scalar field. 

6. Conclusions 

In this paper, we have studied the dynamics of a scalar field in the Lorentz violating scalar- 
vector-tensor theory of gravity, taking into account the effect of the power-law effective 
coupling vector. Since the effective coupling vector be dynamics variable, the equation of 
state is dependent on the coupling parameter. For the model with the power law Hubble 
parameter and coupling vector, we find an exact solution of the equation of state. A 
constant equation of state corresponds to n = 2 while for n > 2 leads to dynamics equation 
of states. In this case, the scalar fields are completely associated with the Lorentz violation. 

Also, the different form in the coupling vector and the potential models lead to the 
different qualitative evolution in two regimes of inflation. The results show that, for the 
inverse power-law potential, the Hubble parameter decreases during the Lorentz violation 
stage and increases in the standard slow roll stage. For the power-law potential, the Hubble 
parameter increases during the Lorentz violation stage but it decreases in the standard slow 
roll stage. 

From the qualitative study of the dynamical system, we have demonstrated the at- 
tractor behavior of inflation driven by a scalar field in the context of scalar-vector-tensor 
theory of gravity. We have found that there exists the Lorentz violating kinetic domi- 
nated solution and the Lorentz violating potential-kinetic dominated solution, depending 
on the region of the coupling parameter in the simplest Lorentz violating chaotic infla- 
tion model. The quadratic coupling vector and the chaotic potential correspond to the 
constants Ai = A2 = — 2y / m and T± = T2 = 1/2. There are two important results of 
this study, which are different from the scalar-tensor theory of gravity: the condition for 
the accelerating universe, Eq. ( |5.14 ) and the slope, p, Eq. ( 5.16 ). The first one yields 



A1X0 > -y/3/8(cj0 + 1/3). The analysis of the critical points show that we may obtain an 
accelerated expansion provided that the solutions are approaching the Lorentz violation 
kinetic dominated solution with m > 1/6 and approaching the Lorentz violation potential- 
kinetic dominated solution with m < 3/8. When the accelerating condition is satisfied, the 
slope p characterizes the properties of the inflating universe: power-law inflation {p > 0), 
de Sitter inflation (p = 0) and superinflation (p = — \p\ < 0). In other cases, if Ai and A2 are 
constants, one finds that the coupling vector is still quadratic in scalar field, (3 ~ c/> 2 , while 
the potential as a function of scalar field <p is given by a power-law potential, V(<p) ~ 27 
and Ti = 1/2, T2 = 1 — I/27, where 7 = A2/A1. Moreover, in order to obtain dynamical 
evolution of the system, we need to solve Eqs. ( |5.6D and (fT?]) together with Eqs. fl5,4|) and 
( |5.5| ). For a realistic model, the effect of an additional component (matter field) would be 
interesting [23]. 
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Finally, we would like to emphasize that there exists an attractor solution in the Lorentz 
violating scalar-vector-tensor theory of gravity. 
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